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The two strands of a DNA molecule with a repetitive sequence can pair into many different 
basepairing patterns. For perfectly periodic sequences, early bulk experiments of Porschke indicate 
the existence of a sliding process, permitting the rapid transition between different relative strand 
positions [Biophys. Chem. 2 (f974) 83]. Here, we use a detailed theoretical model to study the 
basepairing dynamics of periodic and nearly periodic DNA. As suggested by Porschke, DNA sliding 
is mediated by basepairing defects (bulge loops), which can diffuse along the DNA. Moreover, a 
shear force / on opposite ends of the two strands yields a characteristic dynamic response: An 
outward average sliding velocity v ~ f/A is induced in a double strand of length A*', provided / 
is larger than a threshold fc- Conversely, if the strands are initially misaligned, they realign even 
against an external force f < fc- These dynamics effectively result in a viscoelastic behavior of 
DNA under shear forces, with properties that are programmable through the choice of the DNA 
sequence. We Hnd that a small number of mutations in periodic sequences does not prevent DNA 
sliding, but introduces a time delay in the dynamic response. We clarify the mechanism for the 
time delay and describe it quantitatively within a phenomenological model. Based on our hndings, 
we suggest new dynamical roles for DNA in artificial nanoscale devices. The basepairing dynamics 
described here is also relevant for the extension of repetitive sequences inside genomic DNA. 



Introduction 



The basic double-helical structure of DNA is insensi- 
tive to the nucleotide sequence, but many of its biophys- 
ical properties are not. For instance, the local thermo- 
dynamic stability of double stranded DNA (dsDNA) de- 
pends strongly on the sequence Q , and certain sequence 
motifs can cause permanent bends or make DNA more 
bendable Q . Such local modulations of the DNA proper- 
ties play an important role in molecular biology, e.g. for 
nucleosome positioningy] and transcription regulation 
through DNA looping M]. The sequence-dependent sta- 
bility of DNA basepairing is also crucial for applications 
in nanotechnology 0, |E 13 • Clearly, since the thermody- 
namics of DNA basepairing is sequence-dependent, the 
kinetics is sequence-dependent as well. Our aim here is 
to show that the kinetics can display a much richer phe- 
nomenology than might be expected on the basis of the 
thermodynamics alone. 

The dynamics of DNA basepairing can be probed ex- 
perimentally on the single-molecule level with mechanical 
and optical techniques ,9.,. .lalTTlllllTaLll^. One ap- 
proach is to unzip dsDNA from one end of the double 
helix [T^ Il3l [Tsl]. However, unzipping probes only one 
aspect of the basepairing dynamics, the sequential open- 
ing of consecutive basepairs. In a different approach, a 
shear force is applied by grabbing the two strands on 
opposite ends of the dsDNA , see Fig. ^ For a het- 
erogeneous dsDNA with a random sequence, the effect of 
a shear force is to unravel the basepairs from both ends 



jlfil ]. see Fig. which is qualitatively similar to unzip- 
ping. In contrast, with a perfectly periodic sequence, 
e.g. {C)n on the upper and {G)n on the lower strand 
or a higher order repeat such as (CA)jv and {GT)m, the 
two strands can bind in many conhgurations 17]. An 
applied shear force then facilitates local strand slippage 
and can induce macroscopic DNA shding [13, see Fig.^b. 
[Throughout this paper, we use the term 'DNA slippage' 
for microscopic events where a few bases at the end of the 
double strand unbind and rebind shifted by one or sev- 
eral repeat units. In contrast, 'DNA sliding' refers to an 
average large-scale movement of the two strands against 
each other.] 

DNA slippage is an aspect of DNA basepairing dynam- 
ics, which plays an important role in the generation of a 
class of genetic diseases [isj, |2fl| . If local DNA slippage 
occurs in an Okazaki fragment during DNA replication, 
trinucleotide repeats inside genes can get extended be- 




* Electronic address: 
t Electronic address: 



Richard. Neher@physik.lmu.de| 



Ulrich.GerlandQphysik.lmu.del 



FIG. f: DNA under a shear force, (a) A non-periodic se- 
quence unravels from both ends, driven by the length gain 
from converting stacked bases into longer single strand, (b) 
A periodic DNA sequence can open by sliding, mediated by 
bulge loops that are created at the ends and diffuse freely 
along the DNA. When a bulge loop reaches the opposite end, 
the two strands have effectively slipped against each other by 
a distance equal to the loop size. 
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yond a threshold length for the onset of Huntington's 
and other diseases. Such slippage events are possible only 
within the time window that DNA polymerase needs to 
fill in the Okazaki fragment. Thus, the kinetics of strand 
slippage is an important determinant for the frequency 
of repeat extensions. 

We propose that the dynamics of periodic and nearly 
periodic DNA is interesting also for the design of DNA- 
based nanodevices. Indeed, DNA is becoming increas- 
ingly popular as a building block for the assembly of 
nanoscale structures and devices 0, H, ■ These appli- 
cations already exploit the specificity of the basepairing 
interaction, e.g. to direct the assembly of DNA strands 
into predefined architectures, and the dynamics of DNA 
branch migration, e.g. to replace a bound DNA strand by 
a different strand. In the discussion section, we consider 
several new possible applications of DNA in nanotech- 
nology, based on the dynamic properties identified in the 
main part of this paper. 

Finally, DNA sliding is interesting also from a purely 
theoretical perspective. Since simultaneous slippage of 
all basepairs is kinetically inhibited by an extensive acti- 
vation barrier, the macroscopic sliding of DNA strands is 
a complex process involving the dynamics of many local 
basepairing defects il7„ JJJ . The most likely defects are 
bulge loops, see Fig. which are created at the ends 
of the dsDNA, or, in pairs, anywhere along the molecule. 
Once formed, bulge loops diffuse freely along a periodic 
dsDNA until they annihilate with a loop on the oppo- 
site strand or are absorbed at an end. Mutations in the 
periodicity of the DNA sequence create obstacles for the 
diffusion of bulge loops. Effectively, the bulge loop dy- 
namics can be regarded as a reaction-diffusion process of 
particles and antiparticles in one dimension. DNA shear- 
ing experiments render certain aspects of these dynamics 
observable and permit a quantitative characterization. 

The outline of this paper is as follows. First, we de- 
scribe our theoretical model for the energetics and dy- 
namics of DNA basepairing under a shear force. We 
then show that our model leads to the following predic- 
tions: (i) For periodic dsDNA, the combination of poly- 
mer mechanics and basepairing dynamics gives rise to a 
viscoelastic response to shear forces above a threshold 
/c, where both /c and the viscosity rj are programmable 
over a wide range through the DNA sequence. The vis- 
coelastic behavior can be described with the help of a 
mechanical analog model, (ii) DNA sliding is possible 
even when the exact sequence periodicity is destroyed by 
a few mutations, (iii) The mutations affect the viscoelas- 
tic behavior by introducing a programmable time delay 
before sliding commences after a sudden force jump, (iv) 
The mechanism for the time delay can be understood 
within a phenomcnological model, which also permits a 
quantitative description of the full distribution of time 
delays. Taken together, we find that the sequence depen- 
dence of the basepairing dynamics allows to adjust the 
mechanical response of DNA under a shear force over a 
broad range of behaviors. In the last section, we discuss 



the experimental ramifications of these findings. 



DNA model 

To study the dynamics of DNA sliding, we consider 
a DNA molecule under a shear force /, which can be 
applied experimentally by pulling the opposite 5' ends 
16] or, alternatively, the opposite 3' ends. In a coarse- 
grained description, the configuration of the DNA is spec- 
ified by its basepairing pattern S and the spatial con- 
tours of both strands. A generic configuration consists 
of two unstretched and two stretched single strands, and 
the central region from the first to the last basepair, see 
Fig. [2 

We will not explicitly describe the dynamics of the spa- 
tial polymer degrees of freedom, but assume rapid equili- 
bration compared to the timescale of DNA sliding. This 
assumption is justified for short DNA molecules: The 
timescale to equilibrate a semiflexible polymer of length 
L and persistence length Ip in a solvent of viscosity rj is 
vL^/72]Jf, where / is an external force applied to its 
ends |2lj . For a DNA of 150 bp (one persistence length) 
in water at a 10 pN load, the equilibration time is on the 
order of 0.01 /is, which is fast compared to the millisecond 
timescale of DNA sliding observed in the reannealing ex- 
periments of Porschke [13. Hence, we integrate out the 
contour conformations to obtain a free energy function 
E{S) that depends only on the basepairing pattern S. 
The total free energy E{S) can be split up into three 
terms, 

E{S) — Estretch{S) + E\yp{S) + E\oop{S) , (1) 

corresponding to the stretching energy, the free energy 
gain due to basepairing, and the free energy cost of (in- 
ternal or bulge) loops in the pattern 5, respectively. 

Polymer model. — The mechanical polymer properties 
of DNA enter only into the stretching energy, which we 
write in the form Estretch{S) = —fL(S), with an effec- 
tive force-dependent total length L{S) for the stretched 
DNA, i.e. the two single-stranded ends where the force 
is applied and the central DNA segment from the first 
to the last base pair, see Fig. 13 The unstretched sin- 
gle strands do not contribute to the free energy, since 
we take all energies relative to the unstretched and un- 
paired state, which is the usual convention 0. For the 
stretched single strands, we use a freely jointed chain 
polymer model with a Kuhn length twice the bare seg- 
ment length hs ~ 0.7 nm for a single base 0|. With this 
model, each unbound base at the ends where the force is 
applied contributes an effective length bs{f) to the total 
length L{S), where 
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and ksT is the thermal energy. Note that bs{f) differs 
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FIG. 2: DNA free energy model. The free energy E{S) of a 
basepairing pattern S contains three separate contributions: 
(i) A negative binding energy for basepairing. For simplicity 
we assign the same binding energy Eb''*' for every basepair of 
type k, regardless of the neighboring bases, (ii) A positive 
free energy cost for internal and bulge loops. We assign the 
same cost ei for every loop, regardless of its length and base 
sequence, since the detailed choice of the loop cost function 
does not affect our main findings, (iii) A stretching energy. 
For a given pattern S, the stretching energy can be written 
in the form —fL{S), with an effective length L{S), which is 
obtained from force-dependent base-to-base distances 5d(/) 
and bs{f) for double and single strand, respectively. Note 
that L{S) does not correspond to the physical length of the 
DNA molecule (see main text). 



from the average extension of one segment in the direc- 
tion of the force. Instead, the average total extension (x) 
of a DNA with basepairing pattern S is calculated as the 
force derivative of the stretching free energy, 

dEstretch{S) 

i^) = qJ ' (3) 

which yields the correct (Langevin) form for the exten- 
sion of a freely jointed chain. For the central DNA re- 
gion from the first to the last basepair, we assume a B- 
DNA conformation and use a wormlike chain model with 
persistence length Ip = 50 nm and a contour length of 
6d = 0.34 nm per base. [The length of an asymmetric 
loop in the central region is approximated by counting 
only the bases in the shorter arm of the loop.] For the 
forces of interest here, the effective length of a basepair, 
6d(/), is given by the asymptotic formula 

The force-dependence of the lengths bs{f), ^d(/) is in 
fact essential only for our calculation of the viscoelastic 
response. For all other properties considered below, the 
force-dependence has no qualitative effect, and will hence 
be neglected (i.e. bs{f) = bs, bd{f) = &d everywhere 
except in the section 'Viscoelastic behavior'). 

Base-pairing energy model. — To obtain a compact the- 
oretical description, we use a base-pairing energy model 
which is simplified from the nearest-neighbor model of 
Ref. 0, but nevertheless permits semi-quantitative pre- 
dictions. We exclude base pairs within a strand, and 
assign a binding (free) energy Sb > for each base 
pair of type k (Watson- Crick or other) between strands. 



see Fig. 12 Hence, Ehp{S) = - X^fc '^'=('5) ^b^'^\ where 
nj; (S) is the total number of type k basepairs in the base- 
pairing pattern S. Similarly, we assign a loop initiation 
cost Ei > for each internal or bulge loop in a given 
pattern (we neglect an additional length-dependent loop 
cost, which has no qualitative effect on the results dis- 
cussed below). Therefore, E\oop{S) = q{S) ee, where q{S) 
is the total number of loops in the pattern. The numeri- 
cal values of the free energy parameters eb^*'-' and e£ are 
temperature dependent jj]. The actual values used in 
our simulations are given below. 

Elementary kinetic steps. — We support our phe- 
nomenological theory presented below by simulating the 
DNA base-pairing dynamics in detail. To this end, we 
use a kinetic Monte Carlo scheme with three single base 
moves as elementary steps |2^: basepair opening, base- 
pair closing, and 'basepair slippage'. Here, basepair slip- 
page refers to a local shift of the binding partner of a 
base, which is possible only for basepairs next to un- 
bound bases, i.e. inside loops or at the ends of the 
molecule. Clearly, basepair slippage can also be gener- 
ated by a basepair opening move followed by a basepair 
closing move. However, the work of Porschke 0| in- 
dicates that basepair slippage is faster than would be 
expected from the individual rates for basepair opening 
and closing, see below. Hence, we include the basepair 
slippage move, as has been done previously [23| . 

Kinetic rates. — To fully specify our model for the 
DNA base-pairing dynamics, we need to assign a rate 
to each elementary kinetic step. Careful relaxation ex- 
periments determined the rate for basepair closing at 
the end of helical segments to be 1 — 20 x 10^ s~^, where 
the range indicates the experimental uncertainty. In our 
model, we assume that this rate is independent of the 
identity of the basepair. To reproduce the correct equilib- 
rium behavior from our base-pairing dynamics, the rate 
for opening a basepair of type k at the end of a helix 
must be reduced by a factor exp(— eb'^^V^s^) ^i^^ re- 
spect to the closing rate. From reannealing experiments 
with periodic sequences, Porschke estimated the rate 
for the displacement of a bulge loop by one base, i.e. 
the rate for basepair shppage, to be roughly 5 x 10^ s~^. 
Hence, the rates for basepair closing and basepair slip- 
page are roughly equal, within experimental accuracy. In 
our model, we set them exactly equal, for simplicity (our 
main results are in fact insensitive to the precise value of 
the closing rate, see below). In general terms, our model 
assumes that all kinetic rates of passing from a higher 
energy configuration to one with lower energy are the 
same, while the reverse rates are chosen to obey detailed 
balance. It may be noteworthy that recent theoretical 
work on the kinetics of force-induced RNA unfolding, 
which used similar assumptions, produced surprisingly 
good agreement with experiment j25| . 

Below, we report all of our kinetic simulation data in 
units of Monte Carlo steps (MCS). From Ref. [l3|, our 
best estimate for the real time equivalent of one Monte 
Carlo step is 0.2 /zs. However, it should be kept in mind 
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that this estimate bares a large uncertainty. 



Sliding dynamics of periodic sequences 

The simplest periodic sequence is a repetition of one 
base on one strand, e.g. AAA. . ., and the complemen- 
tary base on the other. In this case, we have only base- 
pairs of one type (i.e. Eb^''' = Eb in our model) and each 
base on one strand can bind to any base on the other 
strand. For longer repeat units, e.g. triplet repeats such 
as CAGCAG . . . , which play an important role in genetic 
diseases 19], one can treat a repeat unit as an effective 
base with larger associated energies e^,, eg and lengths 
bd, bg. We are interested in the basepairing dynamics in- 
duced by a constant shear force / that is suddenly turned 
on at i = 0. In the following, we first review the physi- 
cal description of DNA sliding dynamics, which we have 
established already in Ref. [l^]- We then construct a 
mechanical analog model to characterize the viscoelastic 
response of periodic DNA and its sequence-dependence. 

Quantitative phenomenological description. — As illus- 
trated in Fig. sliding of periodic dsDNA is mediated 
by the creation, diffusion, and annihilation of bulge loops. 
When a force is applied, the diffusion of bulge loops 
within the dsDNA remains unbiased, assuming the force 
does not deform the dsDNA structure significantly (this 
assumption clearly breaks down for forces above the B-S 
transition around 65 pN). In contrast, the force strongly 
affects the rates at which bulge loops are created at the 
ends. When the two DNA strands are misaligned, these 
creation rates are imbalanced, since a bulge at an over- 
hanging end does not reduce the number of basepairs in 
the structure (whereas it does on the opposite end) . This 
imbalance produces a restoring force fc, which can be ob- 
tained approximately by balancing the energy cost 
of breaking a basepair with the gain in stretching energy, 
£h = f ■ (2&s(/) — bd{f))- The restoring force creates an 
average inward drift that realigns the two strands. To 
obtain an outward drift velocity v, i.e. macroscopic slid- 
ing, one needs to overcome the restoring force fc, so that 
u ~ (/ — fc) in the vicinity of fc- Indeed, fc becomes a 
critical force in the thermodynamic sense when the limit 
of a large strand length N is taken and the state where 
the strands are completely separated is excluded. 

At the critical force the rates at which bulge loops are 
produced, are equally large on the overhanging stretched 
and the unstretched ends. The average sliding velocity 
V vanishes, however the bulge loop dynamics still leads 
to a macroscopic diffusion of the two strands relative to 
each other, with a diffusion coefficient D. Interestingly, 
this diffusion coefhcient scales with the total number of 
bases as -D '--^ 1/-^, so that the rupture time r required 
to separate the two strands completely scales as r ^ 
instead of the usual r ~ iV^ for diffusion of a particle over 
a distance N. This scaling of D is due to the fact that 
loops are generated at the ends with a constant rate, but 
only result in a global shift between the strands, if they 



diffuse over a distance iV, either to annihilate at the 
other end or with a loop on the opposite strand. In both 
cases, the probability for an event scales as 1/N. The 
D ~ 1/N scaling occurs also in the reptation problem 
of polymer physics, and indeed the microscopic origin is 
closely related, as motion is mediated by defect diffusion 

Since the production of a loop on the stretched ends 
shortens the molecule, the corresponding production rate 
decreases with /. Hence, the rates of events extending or 
shortening the double stranded region, that are equal for 
f — fc, differ at other forces resulting in a drift. Each of 
these rates, and consequently also the sliding velocity v 
is proportional to 1/iV. From the Einstein relation one 
expects w ~ (/ — fc)D ^ 1/-/V, in agreement with this 
result. With the negative (inward) drift velocity below 
fc, rupture events are driven by rare fluctuations, and 
the rupture time t grows exponentially with N, as is 
characteristic for thermally activated transport over an 
extensive energy barrier. On the other hand, for forces 
larger than fc, the scaling of v leads to rupture times 
increasing as t iV^. 

This scaling holds up to a force /*, above which the 
rupture times grow only linearly with N, due to a dynam- 
ical transition in the opening mode from sliding to unrav- 
eling (i.e. the opening mode of heterogeneous dsDNA). 
For / > /*, it is energetically favorable to break base- 
pairs consecutively from both ends and both strands are 
after breaking N basepairs. Within our model, /* is well 
approximated by the solution of / = eb/[6s(/) — ^d(/)], 
i.e. the balance between the basepairing energy and the 
stretching energy gained by extending the molecule by 
the difference between the length of an unbound base 
and a basepair. [For DNA sequences where this force is 
large enough to deform the DNA structure, in particu- 
lar for /* above the ^ 65 pN of the B-S transition, the 
unraveling mode may not exist.] 

Viscoelastic behavior. — DNA sliding can be regarded 
as a viscous flow of the two strands relative to each other. 
According to the physical picture reviewed above, this 
flow has interesting nonlinear and sequence-dependent 
properties. Since the shear force elicits also an elastic 
response (due to the entropic elasticity of DNA), the be- 
havior of periodic dsDNA is reminiscent of a viscoelastic 
material. Such materials combine solid-like and fluid-like 
mechanical properties when probed by external stress. In 
the following we examine this analogy more closely. 

The mechanical behavior of a typical viscoelastic ma- 
terial can be described by a Zener model ^3]i which is 
constructed e.g. by connecting a Kelvin element (a dash- 
pot in parallel with a spring) in series with a spring, see 
Fig. 121 The Zener model reproduces the two prominent 
characteristics of viscoelastic materials: (i) In a 'creep 
experiment', where a constant stress is suddenly applied, 
an instantaneous elastic strain is followed by a gradual 
creep towards a new equilibrium, (ii) When the strain is 
suddenly increased, the stress rises sharply and then re- 
laxes gradually to an equilibrium value. On a qualitative 
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FIG. 3: Viscoelastic response of periodic DNA. (a) The 
shear force on an 80 bp dsDNA (with two 20 bp ssDNA hnk- 
ers) is switched periodically between /mm = 11 A pN and 
fmax = 19 pN (upper panel). The center and bottom panels 
show the extension-time trace for heterogenous and periodic 
DNA, respectively (energy parameters: ~ 1.11 fc^T and 
ee = 2.8 ksT, roughly corresponding to AT basepairs at 50°C 
The time units are Monte Carlo steps (MCS), the real 
time equivalent of which is discussed in section 'Kinetic rates'. 
The heterogeneous DNA responds only elastically to the force 
jumps, mostly due to stretching of the linkers. The length of 
the periodic DNA shows a similar elastic strain, but in addi- 
tion, the molecule elongates at a finite speed due to sliding, 
since fmax > fc = 16.3 pN. When the molecule is relaxed, we 
find an elastic response and inward sliding, since fmin < fc- 
The length of the periodic DNA fluctuates strongly due to 
loop formation and annihilation, (b) The viscoelastic behav- 
ior displayed by a periodic DNA molecule can be described 
by a generalized Zener model, where harmonic springs are 
substituted by anharmonic elastic elements describing poly- 
mer elasticity and the restoring force fc- The ideal dashpot 
(with viscosity rj) creates the viscous behavior of periodic ds- 
DNA. (c) The response of the above idealized model to the 
same periodic force resembles the average response of periodic 
DNA. 



contraction beyond maximal overlap. In contrast to the 
standard Zener model with harmonic springs, the stress 
in response to a strain, will relax to the value /c, indepen- 
dent of the displacement (within a certain range). Fig.|3Ji 
shows extension-time traces obtained from our model, 
both for a periodic (bottom panel) and a heterogeneous 
DNA (center panel), see caption for parameters. Here, we 
have considered a creep test situation where the force is 
switched periodically between fmin < fc and fmax > fc 
(top panel). Fig.|3|; shows the corresponding behavior of 
the generalized Zener model for comparison. We observe 
that the average behavior of the periodic DNA resem- 
bles that of the generalized Zener model, while the het- 
erogeneous DNA shows only elastic behavior. Of course 
the extension also displays strong thermal fluctuations, 
which play an important role in single molecule dynam- 
ics, and ultimately lead to rupture even below the critical 
force [3. 

Programmability- — The viscoelastic behavior de- 
scribed above relies on the basepairing dynamics 
within the DNA molecule, and is manifestly sequence- 
dependent. This fact makes the mechanical behavior of 
dsDNA under shear force programmable, i.e. both the 
force offset fc and the viscosity 77 can be adjusted through 
sequence composition and length of the dsDNA. Even for 
perfectly periodic sequences, there is still a considerable 
freedom in the choice of the sequence composition, since 
a repeat unit can be several bases long and involve dif- 
ferent combinations of Watson- Crick and other basepairs. 
Exploiting this freedom, the range over which the aver- 
age basepairing energy can be programmed is roughly 
0.5 — AUbT which translates into an equally wide 
range of force offsets|34] fc = £b/(2 6s — 6d)- 



level, periodic dsDNA displays these same characteristics 
in its average behavior: (i) Upon sudden application of 
a constant force / in the range /c < / < /*, the DNA 
rapidly stretches against its entropic elasticity and slowly 
creeps with a viscosity rj that is proportional to the num- 
ber of bases in the double strand. However, it will not 
approach a new equilibrium but rupture eventually, (ii) 
When the extension of the DNA is suddenly increased, 
the tension rapidly rises and then slowly relaxes to the 
critical value fc (provided the initial rise was above fc)- 

The viscoelastic behavior of periodic DNA can be de- 
scribed by a nonlinear generalization of the Zener model, 
see Fig.|3lD, where the Kelvin element effectively describes 
the basepairing dynamics, while the outer elastic element 
accounts for the entropic elasticity of the polymer back- 
bone (consisting of dsDNA, ssDNA, and, if present, link- 
ers to the points of force application). Since the base- 
pairing dynamics of two misaligned complementary pe- 
riodic DNA strands produces a restoring force fc, the 
sliding velocity v is proportional to f — fc- The sliding 
dynamics is thus described by a dashpot in parallel with 
a potential generating the restoring force and preventing 



The velocity of macroscopic strand sliding is deter- 
mined by four factors: (i) the mobility of defects, i.e. the 
rate for bulge loop displacement, (ii) the bulge loop den- 
sity, (iii) the inverse strand length, and (iv) the deviation 
of the force / from the critical force fc- The defect den- 
sity depends sensitively on the base-pairing free energy 
and may vary roughly between 0.001 and 0.2 for different 
repeat lengths and temperatures, leaving great freedom 
to adjust the timescale of DNA sliding. Note that since 
only the bulge loop density and not the individual rates 
for basepair closing and opening influence the sliding ve- 
locity, the rate for bulge loop displacement is the only 
crucial rate parameter in our model. By increasing the 
strand length, the sliding velocity can be made arbitrary 
small, or, equivalently, the viscosity rj can be made ar- 
bitrarily large (77 ^ N). Alternatively, rj is increased 
by using longer repeat units, since 77 grows exponentially 
with the free energy cost of creating a bulge loop. An 
order of magnitude estimate [s^ for the lower bound on 
rj yields ~ 10"'^ pN-s/nm, based on reannealing experi- 
ments with homogeneous oligonucleotides of 10 bps 
With these force- and timescales, DNA sliding should be 
well observable in single-molecule experiments. 
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FIG. 4: Top: Extension-time trace for a perfectly periodic 
DNA oi N = 120 basepairs and the same DNA with 7 mu- 
tations, both under a shear force of / = 12.7 pN (energy 
parameters as in Fig. |^ . Whereas the molecule without mu- 
tations starts sliding almost immediately, the molecule with 
mutations fluctuates about its initial length for some time Tw 
before sliding starts. Bottom: The time trace of the binding 
state (open/closed) for the seven mutated basepairs in the se- 
quence. Each mutated basepair (1 to 7) is unbound wherever 
the corresponding thick horizontal line is broken, and bound 
where the line is shown. Note that the mutated basepairs 
do not open/close independently from each other. Instead, 
a mutated basepair opens only once all mutated basepairs to 
the left or right are already open. The black envelope curves 
emphasize the positions of the outermost bound mutation on 
each side. Their dynamics resembles a (biased) random walk. 
Sliding begins when all mutations are open. 



Periodic DNA with weak sequence disorder 

How is the basepairing dynamics affected when a few 
mutations destroy the perfect periodicity? Fig. 0] (top) 
shows two simulated extension-time traces, one for a per- 
fectly periodic sequence and one with M — 7 equidistant 
mutations (DNA parameters: see caption). Here, we as- 
signed the same binding energy to mutated and original 
basepairs, in order to focus on the effects that mutations 
exert on the basepairing dynamics rather than the en- 
ergetics. Furthermore, we assumed that mutated bases 
can only bind to their ground-state binding partners, i.e. 
mutated bases cannot form basepairs with the original 
bases and all mutations are of a different type. The less 
generic effects that can result without these assumptions 
are discussed below. 

We observe that the mutations have a drastic effect: 
whereas the original sequence begins to slide almost im- 
mediately after application of the force, the mutated se- 
quence exhibits a pronounced delay before sliding sets in. 



Indeed, the figure suggests that the mutated sequence has 
two characteristic timescales, a waiting time Tw during 
which the extension fluctuates around a constant value, 
and a sliding time Tg during which the extension increases 
until the two strands are completely separated. A sec- 
ond, less drastic effect of the mutations is to reduce Tg, i.e. 
once sliding starts, it is faster than without mutations. 
Note that the convex shape of both sliding curves is due 
to the fact that the sliding velocity increases as the length 
of the double stranded region decreases, v ^ r]~^ ~ N~^, 
see above. 

What is the physical mechanism that sets the waiting 
timescale ? Clearly, sliding can begin only after all mu- 
tated basepairs have been broken, since otherwise the two 
strands are locked into one relative position. Arguably 
the simplest scenario would be that all mutations inde- 
pendently fluctuate between the open and closed state, 
and sliding commences when all mutations happen to 
be open simultaneously. Alternatively, the dynamics of 
the mutated basepairs could be correlated. To clarify 
the dynamical mechanism, we plotted the binding state 
(bound/unbound) of all mutated basepairs alongside the 
trajectories in Fig. ^(bottom, gray curves). It is evident 
that the mutations do not open independently. Instead, 
interior mutations open only once the neighboring muta- 
tion towards the exterior has already opened. 

Two random walker model. — Inspection of Fig.2|(bot- 
tom) suggests that the positions of the two outermost 
bound mutations might in fact perform a (biased) ran- 
dom walk, see the black curves in Fig. 21 If true, the 
waiting time could be interpreted as the first colli- 
sion time Tc of two random walkers (2RW) on a row of 
M -\- I discrete sites, with force-dependent in- and out- 
ward hopping rates kin and fcout- To test this hypothesis, 
we compare the histogram of r^'s (from many simula- 
tions) with the distribution P(tc) of first collision times 
for 2RW. Fig. [SI shows three such histograms (main panel 
and two insets) obtained with the same DNA parameters 
as in Fig. 13 but with three different forces. Superimposed 
are the distributions P{tc), calculated as described below 
and in Supporting Material. The case shown in the main 
panel of Fig. |31 corresponds to a force value for which the 
2RW are unbiased, i.e. fcjn = fcout, whereas the left inset 
corresponds to a smaller force producing a bias to the 
outside (fcout > fcin) and the right inset shows the oppo- 
site case of a larger force and fcout < fcin- In all three 
cases, the observed histogram is well described by the 
2RW model. Indeed, despite some caveats (see below), 
this model can serve as a useful coarse-grained descrip- 
tion for the basepairing dynamics preceding the sliding 
stage. 

The calculation of the first collision time distribution 
P{tc) belongs to the class of first passage problems, which 
has been studied extensively in statistical physics [28l| . 
In the context of the helix-coil transition, Schwarz and 
Poland 29] (see also |33|) already solved the associated 
diffusion problem. Here, we use their work as a basis 
to treat the first passage problem. One can replace the 
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FIG. 5: Waiting time distributions. Main panel: The his- 
togram of waiting times r„ of a 120 bp long DNA sequence 
with M = 7 equidistant mutations subject to a force / — 
/c — 12.9 pN, is well described by the distribution of collision 
times (dashed line) of the two random walker (RW) model, 
see main text and Fig. |S| The solid line shows the parameter- 
free asymptotic distribution of Eq. Q for comparison. Insets: 
Distribution of Tw for forces above and below /c (/ = 15.2 pN 
and / — 11.4 pN, respectively). The dashed lines are fits 
using the RW model with directional bias, see main text. 



problem of 2RW in one dimension by the equivalent prob- 
lem of one RW on a two dimensional lattice with a tri- 
angular shape, see Fig. El In the following, the unbiased 
case (fcin — fcout = k) is of particular interest. In this 
case, there is only the single rate constant k, which can 
be absorbed in the unit of time, so that the distribution 
P{tc) depends only on the number of lattice points (i.e. 
the number of mutations). However, in the limit of large 
M this dependence also disappears, if we use the rescaled 
collision time fc = Tcfc/M^. The resulting parameter- free 
distribution can be expressed in the form (see Supporting 
Material) 



P{rc) 
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where Q{t) is the rapidly converging series 
^ (-1) 
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This distribution is plotted as the solid line in the main 
panel of Fig. [3 Even when M is small, the distribution 
(O is a good approximation to the actual distribution, as 
illustrated by the dashed line in Fig. O showing the exact 
distribution for the case of M = 7. 

In the case of biased RW (fcin ^ kout), we compute 
P{tc) numerically. The dashed curves in the insets of 
Fig. El show these distributions for M = 7 mutations, 
where we have used the rates fcin, fcout as fit parameters. 

Scaling of mean waiting time. — In the 2RW model. 
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FIG. 6: On a coarse-grained level, the dynamics of mutation 
opening/closing can be described by a model of two parti- 
cles hopping on a ID lattice, with inward/outward hopping 
rates fcin, fcout- Their positions represent the two outermost 
closed mutations. When the particles collide, all mutations 
have opened. Equivalently, one can consider a single particle 
hopping on a triangular 2D lattice. The first collision time 
then corresponds to the time to reach the diagonal absorbing 
boundary. 



the mean first collision time follows the diffusive behav- 
ior Tc ~ when the RW are unbiased, see above. When 
the walkers have an inward bias, this changes to lin- 
ear scaling Tc ~ M, whereas Tc increases exponentially 
with M for an outward bias, see Supporting Material. To 
test these predictions of the 2RW model, we determined 
the mean waiting time (tw) for different M and different 
forces / from our DNA simulations. Fig. [7^ shows (t„) 
as a function of AI (on a double logarithmic scale) for the 
same three force values as in Fig. |S1 Here, we increased 
the total DNA length N proportional to M, in order to 
keep the mutation density constant and equal to that of 
Fig. [SI At the smallest force, the waiting time increases 
exponentially with M, as expected. At the intermediate 
force, corresponding to the case of unbiased walkers, we 
find a scaling (r^) ~ with ( ss 2.4, whereas C ~ 1-5 
for the largest force. We expect that the values of these 
exponents are strongly influenced by finite size effects, 
since we can vary M only over roughly one decade. How- 
ever, our results indicate that the waiting times increase 
more rapidly with the system size than expected on the 
grounds of our phenomenological 2RW model. A possible 
explanation is given in 'Microscopic mechanism' below. 

How does the mean waiting time depend on the ap- 
plied force? Fig.jTlj shows three curves of (tw) vs. / for 
different mutation densities. The vertical dashed lines 
indicate the force value where fcin = fcout for each curve. 
Below these values, (r^) increases sharply with decreas- 
ing force. Indeed, it is reasonable to consider the force 
fc where fcin = fcout as a generalization of the critical 
force fc to the case of weakly disordered sequences. As 
explained in Supporting Material, the rates fcin, fcout can 
be extracted in several different ways from the simulation 
data, leading to fc values which are mutually compatible. 

Fig-Cfc summarizes the different dynamical regimes as 
a function of the applied force / and the mutation den- 
sity v. Without mutations (z/ — 0) the force axis is di- 
vided into three regimes, with rupture driven by rare 
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fluctuations, continuous sliding, and unravelling at low, 
intermediate, and large forces, respectively. As muta- 
tions are introduced (i^ > 0), the boundary fc{i^) between 
the fluctuation-driven 'Kramers regime' and the 'sliding 
regime' rises to larger forces, and the sliding regime ac- 
quires the time delay of Fig. 01 as a new feature. It is 
clear from Fig. [3: that the force interval displaying slid- 
ing behavior becomes narrower as the mutation density 
is increased. This trend can be understood within a more 
microscopic picture, see below. We could not determine 
unambiguously whether the sliding regime vanishes com- 
pletely already at a finite mutation density. However, 
it is clear that sliding will in practice be unobservable 
for sequences with many mutations. The qualitative fea- 
tures depicted in Fig. [3; are robust against variations 
in our microscopic parameters eb, £i, bs{f), and fcd(/)- 
However, the positions of the boundaries between the 
different regimes depend on these parameters, see below. 

Microscopic mechanism. — Why does the mutation dy- 
namics of Fig. ^ (bottom) resemble the behavior of two 
random walkers? First, the opening of a mutation (and 
subsequent local shift of the two strands against each 
other) is always associated with the formation of two per- 
manent loops, see Fig.EJleft). Hence, the opening of mu- 
tations is energetically expensive and mutations remain 
mostly closed, as long as this cost is not compensated 
by any gain in entropy or stretching energy. Since there 
is no such gain when an interior mutation opens, mu- 
tations can only open beginning from the ends towards 
the inside: loops are constantly created at the ends of 
both strands and propagate inwards until they hit a mu- 
tation, which forms a barrier to bulge loop diffusion. On 
the unstretched strand, loops are generated at a higher 
rate than on the stretched strand, resulting in a larger 
quasi-equilibrium loop density 'l8'| . When the outermost 
bound mutation opens spontaneously, the accumulated 
loops on the exterior unstretched strand can suddenly 
penetrate to the inside. This penetration results in an 
entropy gain and a relative shift of the mutated bases, 
which prevents immediate recombination, see Fig. 11 in 
Supporting Material. The size of the entropy gain and 
the shift increases with the distance to the next muta- 
tion. Therefore, the mutation density, not the absolute 
number, is the relevant parameter that determines the 
relative magnitude of the hopping rates fcjn and fcout in 
the random walker model, and hence fixes the value of 
the force threshold /c. 

We now discuss how parameter changes affect the lo- 
cation of the boundaries between the different dynamical 
regimes in Fig. Et. First, it is clear that increasing the 
basepairing energy eb, will shift both force thresholds, 
fc and /*, towards higher forces. Furthermore, from the 
above microscopic picture, it follows that the inward hop- 
ping rate km is proportional to the average loop density, 
while the outward hopping rate /cout decreases with the 
loop density. The average loop density in turn is af- 
fected by our energy parameters: with increasing Eb, ££, 
the average loop density decreases, and consequently the 
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FIG. 7: (a) The mean waiting time (tw) as a function of the 
system size (the mutation density of = 1/15 is kept fixed 
as the number of evenly spaced mutations M is increased). 
At low forces the scaling is exponential (circles: data for / = 
11.4 pN; solid line: exponential fit), while we find power law 
behavior at the force threshold (/c = 12.9 pN, squares) and 
above (/ = 15.2 pN, diamonds), (b) The mean waiting time 
as a function of the applied force for a sequence of A' = 240 bp 
with 5, 9 and 15 mutations. The dashed lines indicate the 
threshold force /c for each case. Below the threshold, (tw) 
rises sharply, (c) Different regimes of the DNA dynamics 
in the parameter space {f,i^). The Kramers regime (DNA 
rupture becomes exponentially slow with increasing system 
size) is separated from the (delayed) sliding regime by the 
line /c(j^) where the inward and outward hopping rates are 
equal, ki„ — kout (circles: data, solid line: interpolation). At 
forces larger than /*, the molecule dissociates by unravelling 
from both ends. 



boundary /c(j^) is shifted towards lower mutation den- 
sities, i.e. the sliding regime becomes more sensitive to 
mutations (this tendency is enhanced by the rising ener- 
getic cost for opening a mutation). 

So far, we considered only mutations with binding en- 
ergy equal to the original bases. Dropping this restriction 
leads to a sloped boundary f*{i^) between the sliding and 
unravelling regimes, and also affects the slope of fc{i^). 
Furthermore, we assumed above that all mutations are 
of a different type and bind only to their native binding 
partner. Without this assumption, bases belonging to 
different mutated basepairs can bind on encounter dur- 
ing the sliding phase. These basepairs have to be opened 
in the same way as during the waiting phase preceding 
sliding. When mutations are equidistant, this effect be- 
comes particularly strong, leading to additional intervals 
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of constant length, i.e. plateaus in the extension vs. time 
trace. Another important effect, caused both by variable 
spacing and energies of mutations, is that the hopping 
rates fcin and fcout become site dependent, so that the 
random walks are effectively on a rugged energy land- 
scape [sills^ . 

Finally, we stress that the 2RW model is phenomeno- 
logical and fails to describe certain features of the DNA 
dynamics. For instance, our simple description has ne- 
glected correlations between subsequent hopping steps of 
a RW, see Supporting Material. Short-range correlations 
do not affect the long-time behavior, which may explain 
why our model describes the shape of the waiting time 
distribution accurately, see Fig. |S1 A more drastic ap- 
proximation is that the 2RW model does not account for 
the time required to bring in new loops from the ends to a 
mutation deep inside the dsDNA. The fact that this time 
increases with the length of the DNA may be the cause 
for the waiting time to rise more rapidly with the system 
size than expected from the 2RW model, see Fig. [7| 

Conclusions and Outlook 

The basepairing dynamics in DNA and RNA molecules 
is only beginning to be explored. Here, we have shown 
that even the seemingly simple case of periodic DNA se- 
quences displays rich behavior, which can be revealed by 
applying a shear force. Our main finding is that the mi- 
croscopic dynamics of bulge loop defects endows DNA 
with viscoelastic properties, which can be programmed 
into the sequence. Weak sequence disorder does not 
abolish these properties, but (i) introduces a delay, since 
all mutations have to be broken before DNA sliding be- 
gins, and (ii) effectively narrows the viscoelastic force 
regime. The dynamics of mutation breaking is an in- 
teresting statistical process, with main features that can 
be understood by considering a first passage problem of 
two random walkers. Our theoretical study has led to 
several experimental ramifications. For instance, we pre- 
dict that periodic or nearly periodic DNA responds to 
sudden stress by slowly relaxing its tension to a thresh- 
old value independent of the initial stress (provided the 
DNA is not too short). This stress relaxation process 



cannot occur for heterogeneous DNA. Furthermore, we 
predict that the relaxation velocity is inversely propor- 
tional to the DNA length, so that the timescale of the 
dynamics can be easily adjusted into the range of interest 
for a given experimental setup. We expect the existence 
of the different dynamical regimes shown in Fig. [T]: to 
be independent of our detailed model assumptions. As 
DNA slippage is directly linked to the production rate 
and mobility of bulge loops, single molecule experiments 
on DNA sliding would test our basic understanding of 
basepairing dynamics in DNA. 

The same properties, that make DNA uniquely suited 
for reliably storing genetic information while keeping it 
accessible, permit many applications in nanotechnology 
5]. For instance, dsDNA has been used as reversible 
crosslinker in polymer networks to switch between dif- 
ferent mechanical properties 0, and even DNA-only 
networks with specified topologies can be constructed, 
exploiting the specificity of the basepairing interaction 
0. In other applications, short dsDNA molecules served 
as programmable force sensors |Q] using the sequence- 
dependence of the mechanical rupture force, or DNA- 
based nano-machines were constructed on the basis of 
the DNA branch migration mechanism js^. Our results 
render several new applications for DNA in nanotechnol- 
ogy conceivable. For instance, complementary periodic 
ssDNA's could be used as self-tightening connections in 
nanostructures: once two such strands found each other, 
they will slide to maximize their overlap until the ten- 
sion reaches a value fc- Periodic or nearly periodic DNA 
could also serve as a viscoelastic crosslinker in polymer 
networks, which should lead to different material proper- 
ties from those observed in . Similarly, DNA networks 
could also be endowed with viscoelastic properties, and 
(nearly) periodic DNA might even be useful as a pro- 
grammable reference molecule for kinetic measurements. 
Of course, which of these and other possible applications 
will turn out to be useful in the end is unclear at the 
present stage. However, wc feel that there is a clear po- 
tential that should be explored. 
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